The Larmor frequency shift of a white matter magnetic microstructure model with multiple sources.

Magnetic susceptibility imaging may provide valuable information about chemical composition and microstructural organization of tissue. However, its estimation from the MRI signal phase is particularly difficult as it is sensitive to magnetic tissue properties ranging from the molecular to the macroscopic scale. The MRI Larmor frequency shift measured in white matter (WM) tissue depends on the myelinated axons and other magnetizable sources such as iron-filled ferritin. We have previously derived the Larmor frequency shift arising from a dense medium of cylinders with scalar susceptibility and arbitrary orientation dispersion. Here, we extend our model to include microscopic WM susceptibility anisotropy as well as spherical inclusions with scalar susceptibility to represent subcellular structures, biologically stored iron, and so forth. We validate our analytical results with computer simulations and investigate the feasibility of estimating susceptibility using simple iterative linear least squares without regularization or preconditioning. This is done in a digital brain phantom synthesized from diffusion MRI measurements of an ex vivo mouse brain at ultra-high field.

Rat superior colliculus encodes the transition between static and dynamic vision modes.

The visual continuity illusion involves a shift in visual perception from static to dynamic vision modes when the stimuli arrive at high temporal frequency, and is critical for recognizing objects moving in the environment. However, how this illusion is encoded across the visual pathway remains poorly understood, with disparate frequency thresholds at retinal, cortical, and behavioural levels suggesting the involvement of other brain areas. Here, we employ a multimodal approach encompassing behaviour, whole-brain functional MRI, and electrophysiological measurements, for investigating the encoding of the continuity illusion in rats. Behavioural experiments report a frequency threshold of 18±2 Hz. Functional MRI reveal that superior colliculus signals transition from positive to negative at the behaviourally-driven threshold, unlike thalamic and cortical areas. Electrophysiological recordings indicate that these transitions are underpinned by neural activation/suppression. Lesions in the primary visual cortex reveal this effect to be intrinsic to the superior colliculus (under a cortical gain effect). Our findings highlight the superior colliculus' crucial involvement in encoding temporal frequency shifts, especially the change from static to dynamic vision modes.

Extensive topographic remapping and functional sharpening in the adult rat visual pathway upon first visual experience.

Understanding the dynamics of stability/plasticity balances during adulthood is pivotal for learning, disease, and recovery from injury. However, the brain-wide topography of sensory remapping remains unknown. Here, using a first-of-its-kind setup for delivering patterned visual stimuli in a rodent magnetic resonance imaging (MRI) scanner, coupled with biologically inspired computational models, we noninvasively mapped brain-wide properties-receptive fields (RFs) and spatial frequency (SF) tuning curves-that were insofar only available from invasive electrophysiology or optical imaging. We then tracked the RF dynamics in the chronic visual deprivation model (VDM) of plasticity and found that light exposure progressively promoted a large-scale topographic remapping in adult rats. Upon light exposure, the initially unspecialized visual pathway progressively evidenced sharpened RFs (smaller and more spatially selective) and enhanced SF tuning curves. Our findings reveal that visual experience following VDM reshapes both structure and function of the visual system and shifts the stability/plasticity balance in adults.

MP-PCA denoising of fMRI time-series data can lead to artificial activation "spreading".

MP-PCA denoising has become the method of choice for denoising MRI data since it provides an objective threshold to separate the signal components from unwanted thermal noise components. In rodents, thermal noise in the coils is an important source of noise that can reduce the accuracy of activation mapping in fMRI. Further confounding this problem, vendor data often contains zero-filling and other post-processing steps that may violate MP-PCA assumptions. Here, we develop an approach to denoise vendor data and assess activation "spreading" caused by MP-PCA denoising in rodent task-based fMRI data. Data was obtained from N = 3 mice using conventional multislice and ultrafast fMRI acquisitions (1 s and 50 ms temporal resolution, respectively), using a visual stimulation paradigm. MP-PCA denoising produced SNR gains of 64% and 39%, and Fourier Spectral Amplitude (FSA) increases in BOLD maps of 9% and 7% for multislice and ultrafast data, respectively, when using a small [2 2] denoising window. Larger windows provided higher SNR and FSA gains with increased spatial extent of activation that may or may not represent real activation. Simulations showed that MP-PCA denoising can incur activation "spreading" with increased false positive rate and smoother functional maps due to local "bleeding" of principal components, and that the optimal denoising window for improved specificity of functional mapping, based on Dice score calculations, depends on the data's tSNR and functional CNR. This "spreading" effect applies also to another recently proposed low-rank denoising method (NORDIC), although to a lesser degree. Our results bode well for enhancing spatial and/or temporal resolution in future fMRI work, while taking into account the sensitivity/specificity trade-offs of low-rank denoising methods.

To mask or not to mask? Investigating the impact of accounting for spatial frequency distributions and susceptibility sources on QSM quality.

PURPOSE: Estimating magnetic susceptibility using MRI depends on inverting a forward relationship between the susceptibility and measured Larmor frequency. However, an often-overlooked constraint in susceptibility fitting is that the Larmor frequency is only measured inside the sample, and after successful background field removal, susceptibility sources should only reside inside the same sample. Here, we test the impact of accounting for these constraints in susceptibility fitting. THEORY AND METHODS: Two different digital brain phantoms with scalar susceptibility were examined. We used the MEDI phantom, a simple phantom with no background fields, to examine the effect of the imposed constraints for various levels of SNR. Next, we considered the QSM reconstruction challenge 2.0 phantom with and without background fields. We estimated the parameter accuracy of openly-available QSM algorithms by comparing fitting results to the ground truth. Next, we implemented the mentioned constraints and compared to the standard approach. RESULTS: Including the spatial distribution of frequencies and susceptibility sources decreased the RMS-error compared to standard QSM on both brain phantoms when background fields were absent. When background field removal was unsuccessful, as is presumably the case in most in vivo conditions, it is better to allow sources outside the brain. CONCLUSION: Informing QSM algorithms about the location of susceptibility sources and where Larmor frequency was measured improves susceptibility fitting for realistic SNR levels and efficient background field removal. However, the latter remains the bottleneck of the algorithm. Allowing for external sources regularizes unsuccessful background field removal and is currently the best strategy in vivo.

Temporal Diffusion Ratio (TDR) for imaging restricted diffusion: Optimisation and pre-clinical demonstration.

Temporal Diffusion Ratio (TDR) is a recently proposed dMRI technique (Dell'Acqua et al., proc. ISMRM 2019) which provides contrast between areas with restricted diffusion and areas either without restricted diffusion or with length scales too small for characterisation. Hence, it has a potential for informing on pore sizes, in particular the presence of large axon diameters or other cellular structures. TDR employs the signal from two dMRI acquisitions obtained with the same, large, b-value but with different diffusion gradient waveforms. TDR is advantageous as it employs standard acquisition sequences, does not make any assumptions on the underlying tissue structure and does not require any model fitting, avoiding issues related to model degeneracy. This work for the first time introduces and optimises the TDR method in simulation for a range of different tissues and scanner constraints and validates it in a pre-clinical demonstration. We consider both substrates containing cylinders and spherical structures, representing cell soma in tissue. Our results show that contrasting an acquisition with short gradient duration, short diffusion time and high gradient strength with an acquisition with long gradient duration, long diffusion time and low gradient strength, maximises the TDR contrast for a wide range of pore configurations. Additionally, in the presence of Rician noise, computing TDR from a subset (50% or fewer) of the acquired diffusion gradients rather than the entire shell as proposed originally further improves the contrast. In the last part of the work the results are demonstrated experimentally on rat spinal cord. In line with simulations, the experimental data shows that optimised TDR improves the contrast compared to non-optimised TDR. Furthermore, we find a strong correlation between TDR and histology measurements of axon diameter. In conclusion, we find that TDR has great potential and is a very promising alternative (or potentially complement) to model-based approaches for informing on pore sizes and restricted diffusion in general.

Intrinsic macroscale oscillatory modes driving long range functional connectivity in female rat brains detected by ultrafast fMRI.

Spontaneous fluctuations in functional magnetic resonance imaging (fMRI) signals correlate across distant brain areas, shaping functionally relevant intrinsic networks. However, the generative mechanism of fMRI signal correlations, and in particular the link with locally-detected ultra-slow oscillations, are not fully understood. To investigate this link, we record ultrafast ultrahigh field fMRI signals (9.4 Tesla, temporal resolution = 38 milliseconds) from female rats across three anesthesia conditions. Power at frequencies extending up to 0.3 Hz is detected consistently across rat brains and is modulated by anesthesia level. Principal component analysis reveals a repertoire of modes, in which transient oscillations organize with fixed phase relationships across distinct cortical and subcortical structures. Oscillatory modes are found to vary between conditions, resonating at faster frequencies under medetomidine sedation and reducing both in number, frequency, and duration with the addition of isoflurane. Peaking in power within clear anatomical boundaries, these oscillatory modes point to an emergent systemic property. This work provides additional insight into the origin of oscillations detected in fMRI and the organizing principles underpinning spontaneous long-range functional connectivity.

Larmor frequency shift from magnetized cylinders with arbitrary orientation distribution.

The magnetic susceptibility of tissue can provide valuable information about its chemical composition and microstructural organization. However, the relation between the magnetic microstructure and the measurable Larmor frequency shift is understood only for a few idealized cases. Here we analyze the microstructure formed by magnetized, NMR-invisible infinite cylinders suspended in an NMR-reporting fluid. Through simulations, we scrutinize various geometries of mesoscopic Lorentz cavities and inclusions, and show that the cavity size should be approximately one order of magnitude larger than the width of the inclusions. We also analytically derive the Larmor frequency shift for a population of cylinders with arbitrary orientation dispersion and show that it is determined by the l = 2 Laplace expansion coefficients p 2 m of the cylinders' orientation distribution function. Our work underscores the need to account for microstructural organization when estimating magnetic tissue properties.

Tensor denoising of multidimensional MRI data.

PURPOSE: To develop a denoising strategy leveraging redundancy in high-dimensional data. THEORY AND METHODS: The SNR fundamentally limits the information accessible by MRI. This limitation has been addressed by a host of denoising techniques, recently including the so-called MPPCA: principal component analysis of the signal followed by automated rank estimation, exploiting the Marchenko-Pastur distribution of noise singular values. Operating on matrices comprised of data patches, this popular approach objectively identifies noise components and, ideally, allows noise to be removed without introducing artifacts such as image blurring, or nonlocal averaging. The MPPCA rank estimation, however, relies on a large number of noise singular values relative to the number of signal components to avoid such ill effects. This condition is unlikely to be met when data patches and therefore matrices are small, for example due to spatially varying noise. Here, we introduce tensor MPPCA (tMPPCA) for the purpose of denoising multidimensional data, such as from multicontrast acquisitions. Rather than combining dimensions in matrices, tMPPCA uses each dimension of the multidimensional data's inherent tensor-structure to better characterize noise, and to recursively estimate signal components. RESULTS: Relative to matrix-based MPPCA, tMPPCA requires no additional assumptions, and comparing the two in a numerical phantom and a multi-TE diffusion MRI data set, tMPPCA dramatically improves denoising performance. This is particularly true for small data patches, suggesting that tMPPCA can be especially beneficial in such cases. CONCLUSIONS: The MPPCA denoising technique can be extended to high-dimensional data with improved performance for smaller patch sizes.

Multiple Myeloma-Derived Extracellular Vesicles Modulate the Bone Marrow Immune Microenvironment.

Multiple myeloma (MM), the third most frequent hematological cancer worldwide, is characterized by the proliferation of neoplastic plasma cells in the bone marrow (BM). One of the hallmarks of MM is a permissive BM microenvironment. Increasing evidence suggests that cell-to-cell communication between myeloma and immune cells via tumor cell-derived extracellular vesicles (EV) plays a key role in the pathogenesis of MM. Hence, we aimed to explore BM immune alterations induced by MM-derived EV. For this, we inoculated immunocompetent BALB/cByJ mice with a myeloma cell line, MOPC315.BM, inducing a MM phenotype. Upon tumor establishment, characterization of the BM microenvironment revealed the expression of both activation and suppressive markers by lymphocytes, such as granzyme B and PD-1, respectively. In addition, conditioning of the animals with MOPC315.BM-derived EV, before transplantation of the MOPC315.BM tumor cells, did not anticipate the disease phenotype. However, it induced features of suppression in the BM milieu, such as an increase in PD-1 expression by CD4+ T cells. Overall, our findings reveal the involvement of MOPC315.BM-derived EV protein content as promoters of immune niche remodeling, strengthening the importance of assessing the mechanisms by which MM may impact the immune microenvironment.

Stratification of radiosensitive brain metastases based on an actionable S100A9/RAGE resistance mechanism.

Whole-brain radiotherapy (WBRT) is the treatment backbone for many patients with brain metastasis; however, its efficacy in preventing disease progression and the associated toxicity have questioned the clinical impact of this approach and emphasized the need for alternative treatments. Given the limited therapeutic options available for these patients and the poor understanding of the molecular mechanisms underlying the resistance of metastatic lesions to WBRT, we sought to uncover actionable targets and biomarkers that could help to refine patient selection. Through an unbiased analysis of experimental in vivo models of brain metastasis resistant to WBRT, we identified activation of the S100A9-RAGE-NF-κB-JunB pathway in brain metastases as a potential mediator of resistance in this organ. Targeting this pathway genetically or pharmacologically was sufficient to revert the WBRT resistance and increase therapeutic benefits in vivo at lower doses of radiation. In patients with primary melanoma, lung or breast adenocarcinoma developing brain metastasis, endogenous S100A9 levels in brain lesions correlated with clinical response to WBRT and underscored the potential of S100A9 levels in the blood as a noninvasive biomarker. Collectively, we provide a molecular framework to personalize WBRT and improve its efficacy through combination with a radiosensitizer that balances therapeutic benefit and toxicity.

In vivo Correlation Tensor MRI reveals microscopic kurtosis in the human brain on a clinical 3T scanner.

Diffusion MRI (dMRI) has become one of the most important imaging modalities for noninvasively probing tissue microstructure. Diffusional Kurtosis MRI (DKI) quantifies the degree of non-Gaussian diffusion, which in turn has been shown to increase sensitivity towards, e.g., disease and orientation mapping in neural tissue. However, the specificity of DKI is limited as different sources can contribute to the total intravoxel diffusional kurtosis, including: variance in diffusion tensor magnitudes (Kiso), variance due to diffusion anisotropy (Kaniso), and microscopic kurtosis (µK) related to restricted diffusion, microstructural disorder, and/or exchange. Interestingly, µK is typically ignored in diffusion MRI signal modelling as it is assumed to be negligible in neural tissues. However, recently, Correlation Tensor MRI (CTI) based on Double-Diffusion-Encoding (DDE) was introduced for kurtosis source separation, revealing non negligible µK in preclinical imaging. Here, we implemented CTI for the first time on a clinical 3T scanner and investigated the sources of total kurtosis in healthy subjects. A robust framework for kurtosis source separation in humans is introduced, followed by estimation of µK (and the other kurtosis sources) in the healthy brain. Using this clinical CTI approach, we find that µK significantly contributes to total diffusional kurtosis both in grey and white matter tissue but, as expected, not in the ventricles. The first µK maps of the human brain are presented, revealing that the spatial distribution of µK provides a unique source of contrast, appearing different from isotropic and anisotropic kurtosis counterparts. Moreover, group average templates of these kurtosis sources have been generated for the first time, which corroborated our findings at the underlying individual-level maps. We further show that the common practice of ignoring µK and assuming the multiple Gaussian component approximation for kurtosis source estimation introduces significant bias in the estimation of other kurtosis sources and, perhaps even worse, compromises their interpretation. Finally, a twofold acceleration of CTI is discussed in the context of potential future clinical applications. We conclude that CTI has much potential for future in vivo microstructural characterizations in healthy and pathological tissue.

Short TE downfield magnetic resonance spectroscopy in a mouse model of brain glioma.

PURPOSE: Enhanced cell proliferation in tumors can be associated with altered metabolic profiles and dramatic microenvironmental changes. Downfield magnetic resonance spectroscopy (MRS) has received increasing attention due to its ability to report on labile resonances of molecules not easily detected in upfield 1 H MRS. Image-selected-in-vivo-spectroscopy-relaxation enhanced MRS (iRE-MRS) was recently introduced for acquiring short echo-time (TE) spectra. Here, iRE-MRS was used to investigate in-vivo downfield spectra in glioma-bearing mice. METHODS: Experiments were performed in vivo in an immunocompetent glioma mouse model at 9.4 T using a cryogenic coil. iRE-MRS spectra were acquired in N = 6 glioma-bearing mice (voxel size = 2.23 mm3 ) and N = 6 control mice. Spectra were modeled by a sum of Lorentzian peaks simulating known downfield resonances, and differences between controls and tumors were quantified using relative peak areas. RESULTS: Short TE tumor spectra exhibited large qualitative differences compared to control spectra. Most peaks appeared modulated, with strong attenuation of NAA (∼7.82, 7.86 ppm) and changes in relative peak areas between 6.75 and 8.49 ppm. Peak areas tended to be smaller for DF6.83 , DF7.60 , DF8.18 and NAA; and larger for DF7.95 and DF8.24 . Differences were also detected in signals resonating above 8.5 ppm, assumed to arise from NAD+. CONCLUSIONS: In-vivo downfield 1 H iRE-MRS of mouse glioma revealed differences between controls and tumor bearing mice, including in metabolites which are not easily detectable in the more commonly investigated upfield spectrum. These findings motivate future downfield MRS investigations exploring pH and exchange contributions to these differences.

Soma and Neurite Density MRI (SANDI) of the in-vivo mouse brain and comparison with the Allen Brain Atlas.

Diffusion MRI (dMRI) provides unique insights into the neural tissue milieu by probing interactions between diffusing molecules and tissue microstructure. Most dMRI techniques focus on white matter (WM) tissues, nevertheless, interest in gray matter characterizations is growing. The Soma and Neurite Density MRI (SANDI) methodology harnesses a model incorporating water diffusion in spherical objects (assumed to be associated with cell bodies) and in impermeable "sticks" (assumed to represent neurites), which potentially enables the characterization of cellular and neurite densities. Recognising the importance of rodents in animal models of development, aging, plasticity, and disease, we here employ SANDI for in-vivo preclinical imaging and provide a first validation of the methodology by comparing SANDI metrics with cellular density reflected by the Allen mouse brain atlas. SANDI was implemented on a 9.4T scanner equipped with a cryogenic coil, and in-vivo experiments were carried out on N = 6 mice. Pixelwise, ROI-based, and atlas comparisons were performed, magnitude vs. real-valued analyses were compared, and shorter acquisitions with reduced the number of b-value shells were investigated. Our findings reveal good reproducibility of the SANDI parameters, including the sphere and stick fractions, as well as sphere size (CoV < 7%, 12% and 3%, respectively). Additionally, we find a very good rank correlation between SANDI-driven sphere fraction and Allen mouse brain atlas contrast that represents cellular density. We conclude that SANDI is a viable preclinical MRI technique that can greatly contribute to research on brain tissue microstructure.

Diffusion time dependence, power-law scaling, and exchange in gray matter.

Characterizing neural tissue microstructure is a critical goal for future neuroimaging. Diffusion MRI (dMRI) provides contrasts that reflect diffusing spins' interactions with myriad microstructural features of biological systems. However, the specificity of dMRI remains limited due to the ambiguity of its signals vis-à-vis the underlying microstructure. To improve specificity, biophysical models of white matter (WM) typically express dMRI signals according to the Standard Model (SM) and have more recently in gray matter (GM) taken spherical compartments into account (the SANDI model) in attempts to represent cell soma. The validity of the assumptions underlying these models, however, remains largely undetermined, especially in GM. To validate these assumptions experimentally, observing their unique, functional properties, such as the b-1/2 power-law associated with one-dimensional diffusion, has emerged as a fruitful strategy. The absence of this signature in GM, in turn, has been explained by neurite water exchange, non-linear morphology, and/or by obscuring soma signal contributions. Here, we present diffusion simulations in realistic neurons demonstrating that curvature and branching does not destroy the stick power-law behavior in impermeable neurites, but also that their signal is drowned by the soma signal under typical experimental conditions. Nevertheless, by studying the GM dMRI signal's behavior as a function of diffusion weighting as well as time, we identify an attainable experimental regime in which the neurite signal dominates. Furthermore, we find that exchange-driven time dependence produces a signal behavior opposite to that which would be expected from restricted diffusion, thereby providing a functional signature that disambiguates the two effects. We present data from dMRI experiments in ex vivo rat brain at ultrahigh field of 16.4T and observe a time dependence that is consistent with substantial exchange but also with a GM stick power-law. The first finding suggests significant water exchange between neurites and the extracellular space while the second suggests a small sub-population of impermeable neurites. To quantify these observations, we harness the Kärger exchange model and incorporate the corresponding signal time dependence in the SM and SANDI models.

Glucose fluxes in glycolytic and oxidative pathways detected in vivo by deuterium magnetic resonance spectroscopy reflect proliferation in mouse glioblastoma.

OBJECTIVES: Glioblastoma multiforme (GBM), the most aggressive glial brain tumors, can metabolize glucose through glycolysis and mitochondrial oxidation pathways. While specific dependencies on those pathways are increasingly associated with treatment response, detecting such GBM subtypes in vivo remains elusive. Here, we develop a dynamic glucose-enhanced deuterium spectroscopy (DGE 2H-MRS) approach for differentially assessing glucose turnover rates through glycolysis and mitochondrial oxidation in mouse GBM and explore their association with histologic features of the tumor and its microenvironment. MATERIALS AND METHODS: GL261 and CT2A glioma allografts were induced in immunocompetent mice and scanned in vivo at 9.4 Tesla, harnessing DGE 2H-MRS with volume selection and Marchenko-Pastur PCA (MP-PCA) denoising to achieve high temporal resolution. Each tumor was also classified by histopathologic analysis and assessed for cell proliferation (Ki67 immunostaining), while the respective cell lines underwent in situ extracellular flux analysis to assess mitochondrial function. RESULTS: MP-PCA denoising of in vivo DGE 2H-MRS data significantly improved the time-course detection (~2-fold increased Signal-to-Noise Ratio) and fitting precision (-19 ± 1 % Cramér-Rao Lower Bounds) of 2H-labelled glucose, and glucose-derived glutamate-glutamine (Glx) and lactate pools in GL261 and CT2A orthotopic tumors. Kinetic modeling further indicated inter-tumor heterogeneity of glucose consumption rate for glycolysis and oxidation during a defined epoch of active proliferation in both cohorts (19 ± 1 days post-induction), with consistent volumes (38.3 ± 3.4 mm3) and perfusion properties prior to marked necrosis. Histopathologic analysis of these tumors revealed clear differences in tumor heterogeneity between the two GBM models, aligned with metabolic differences of the respective cell lines monitored in situ. Importantly, glucose oxidation (i.e. Glx synthesis and elimination rates: 0.40 ± 0.08 and 0.12 ± 0.03 mM min-1, respectively) strongly correlated with cell proliferation across the pooled cohorts (R = 0.82, p = 0.001; and R = 0.80, p = 0.002, respectively), regardless of tumor morphologic features or in situ metabolic characteristics of each GBM model. CONCLUSIONS: Our fast DGE 2H-MRS enables the quantification of glucose consumption rates through glycolysis and mitochondrial oxidation in mouse GBM, which is relevant for assessing their modulation in vivo according to tumor microenvironment features such as cell proliferation. This novel application augurs well for non-invasive metabolic characterization of glioma or other cancers with mitochondrial oxidation dependencies.

Correlation Tensor MRI deciphers underlying kurtosis sources in stroke.

Noninvasively detecting and characterizing modulations in cellular scale micro-architecture remains a desideratum for contemporary neuroimaging. Diffusion MRI (dMRI) has become the mainstay methodology for probing microstructure, and, in ischemia, its contrasts have revolutionized stroke management. Diffusion kurtosis imaging (DKI) has been shown to significantly enhance the sensitivity of stroke detection compared to its diffusion tensor imaging (DTI) counterparts. However, the interpretation of DKI remains ambiguous as its contrast may arise from competing kurtosis sources related to the anisotropy of tissue components, diffusivity variance across components, and microscopic kurtosis (e.g., arising from cross-sectional variance, structural disorder, and restriction). Resolving these sources may be fundamental for developing more specific imaging techniques for stroke management, prognosis, and understanding its pathophysiology. In this study, we apply Correlation Tensor MRI (CTI) - a double diffusion encoding (DDE) methodology recently introduced for deciphering kurtosis sources based on the unique information captured in DDE's diffusion correlation tensors - to investigate the underpinnings of kurtosis measurements in acute ischemic lesions. Simulations for the different kurtosis sources revealed specific signatures for cross-sectional variance (representing neurite beading), edema, and cell swelling. Ex vivo CTI experiments at 16.4 T were then performed in an experimental photothrombotic stroke model 3 h post-stroke (N = 10), and successfully separated anisotropic, isotropic, and microscopic non-Gaussian diffusion sources in the ischemic lesions. Each of these kurtosis sources provided unique contrasts in the stroked area. Particularly, microscopic kurtosis was shown to be a primary "driver" of total kurtosis upon ischemia; its large increases, coupled with decreases in anisotropic kurtosis, are consistent with the expected elevation in cross-sectional variance, likely linked to beading effects in small objects such as neurites. In vivo experiments at 9.4 T at the same time point (3 h post ischemia, N = 5) demonstrated the stability and relevance of the findings and showed that fixation is not a dominant confounder in our findings. In future studies, the different CTI contrasts may be useful to address current limitations of stroke imaging, e.g., penumbra characterization, distinguishing lesion progression form tissue recovery, and elucidating pathophysiological correlates.

Optogenetic activation of striatal D1R and D2R cells differentially engages downstream connected areas beyond the basal ganglia.

The basal ganglia (BG) are a group of subcortical nuclei responsible for motor and executive function. Central to BG function are striatal cells expressing D1 (D1R) and D2 (D2R) dopamine receptors. D1R and D2R cells are considered functional antagonists that facilitate voluntary movements and inhibit competing motor patterns, respectively. However, whether they maintain a uniform function across the striatum and what influence they exert outside the BG is unclear. Here, we address these questions by combining optogenetic activation of D1R and D2R cells in the mouse ventrolateral caudoputamen with fMRI. Striatal D1R/D2R stimulation evokes distinct activity within the BG-thalamocortical network and differentially engages cerebellar and prefrontal regions. Computational modeling of effective connectivity confirms that changes in D1R/D2R output drive functional relationships between these regions. Our results suggest a complex functional organization of striatal D1R/D2R cells and hint toward an interconnected fronto-BG-cerebellar network modulated by striatal D1R and D2R cells.

On the generalizability of diffusion MRI signal representations across acquisition parameters, sequences and tissue types: Chronicles of the MEMENTO challenge.

Diffusion MRI (dMRI) has become an invaluable tool to assess the microstructural organization of brain tissue. Depending on the specific acquisition settings, the dMRI signal encodes specific properties of the underlying diffusion process. In the last two decades, several signal representations have been proposed to fit the dMRI signal and decode such properties. Most methods, however, are tested and developed on a limited amount of data, and their applicability to other acquisition schemes remains unknown. With this work, we aimed to shed light on the generalizability of existing dMRI signal representations to different diffusion encoding parameters and brain tissue types. To this end, we organized a community challenge - named MEMENTO, making available the same datasets for fair comparisons across algorithms and techniques. We considered two state-of-the-art diffusion datasets, including single-diffusion-encoding (SDE) spin-echo data from a human brain with over 3820 unique diffusion weightings (the MASSIVE dataset), and double (oscillating) diffusion encoding data (DDE/DODE) of a mouse brain including over 2520 unique data points. A subset of the data sampled in 5 different voxels was openly distributed, and the challenge participants were asked to predict the remaining part of the data. After one year, eight participant teams submitted a total of 80 signal fits. For each submission, we evaluated the mean squared error, the variance of the prediction error and the Bayesian information criteria. The received submissions predicted either multi-shell SDE data (37%) or DODE data (22%), followed by cartesian SDE data (19%) and DDE (18%). Most submissions predicted the signals measured with SDE remarkably well, with the exception of low and very strong diffusion weightings. The prediction of DDE and DODE data seemed more challenging, likely because none of the submissions explicitly accounted for diffusion time and frequency. Next to the choice of the model, decisions on fit procedure and hyperparameters play a major role in the prediction performance, highlighting the importance of optimizing and reporting such choices. This work is a community effort to highlight strength and limitations of the field at representing dMRI acquired with trending encoding schemes, gaining insights into how different models generalize to different tissue types and fiber configurations over a large range of diffusion encodings.

Evidence for microscopic kurtosis in neural tissue revealed by correlation tensor MRI.

PURPOSE: The impact of microscopic diffusional kurtosis (µK), arising from restricted diffusion and/or structural disorder, remains a controversial issue in contemporary diffusion MRI (dMRI). Recently, correlation tensor imaging (CTI) was introduced to disentangle the sources contributing to diffusional kurtosis, without relying on a-priori multi-gaussian component (MGC) or other microstructural assumptions. Here, we investigated µK in in vivo rat brains and assessed its impact on state-of-the-art methods ignoring µK. THEORY AND METHODS: CTI harnesses double diffusion encoding (DDE) experiments, which were here improved for speed and minimal bias using four different sets of acquisition parameters. The robustness of the improved CTI protocol was assessed via simulations. In vivo CTI acquisitions were performed in healthy rat brains using a 9.4T pre-clinical scanner equipped with a cryogenic coil, and targeted the estimation of µK, anisotropic kurtosis, and isotropic kurtosis. RESULTS: The improved CTI acquisition scheme substantially reduces scan time and importantly, also minimizes higher-order-term biases, thus enabling robust µK estimation, alongside Kaniso and Kiso metrics. Our CTI experiments revealed positive µK both in white and gray matter of the rat brain in vivo; µK is the dominant kurtosis source in healthy gray matter tissue. The non-negligible µK substantially were found to bias prior MGC analyses of Kiso and Kaniso . CONCLUSIONS: Correlation Tensor MRI offers a more accurate and robust characterization of kurtosis sources than its predecessors. µK is non-negligible in vivo in healthy white and gray matter tissues and could be an important biomarker for future studies. Our findings thus have both theoretical and practical implications for future dMRI research.